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Abstract 

We discuss the nonfactorization effects in i? —> J/'h + Xg, which is similar to the 
nonperturbative effect found by Voloshin in the decay B —> y + Xg. The QCD sum rule 
has been used to estimate the hadron matrix elements. We find that the correction 
from this effect is very large and the large discrepancy between the theory and the 
experimental data can be reduced considerably. 
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It is amazing that there is a large discrepancy between the theoretical predictions, which 
is based on the operator product expansion (OPE) with assumption of the factorization, and 
experimental data in the process B —> J/'^ + Xs. The theoretical prediction of the branching 
ratio B —> J/T+X^ is 0.23%, which is only one-third of the experiment data (0.8±0.08)%|jl[|. 
In order to resolve this puzzle, some scenarios have been suggested 0. Besides the possible 
existence of new charmonium states above the DD threshold, nonfactorization effects have 
been widely considered. In the later case, color-octet effect seems to be the most reasonable 
source that may enhance the theoretical calculation. Color-octet mechanism was hrst used 
in calculating cross section of cc production, where the theoretical prediction based on color- 
singlet mechanism is also much smaller than the experimental data. In this mechanism, 
the gluon parton has been taken account so that the color-octet operators responsible for 
b ^ cc + q may give a non-zero contribution to the decay mentioned above. However, the 
most important parameter, the matrix element of color-octet operators, is not calculable so 
far. It needs future determination both from various experiments and theoretical estimations. 

Recently, a nonperturbative correction of order of to R —>• 7 -|- is hrst 

discussed by Voloshin [§ ( a simlar correction is also independently discussed by @] ) which 
arises from the contribution of the gluon-photon penguin graph shown in Fig.l, and then 
discussed by Grant et al , Legeti et a/0 and Buchalla et a/0 . It is obvious that this effect 
also exits in R ^ J/T -|- X^ as nonfactorization effects. Although this correction is not 
large in R ^ 7 -|- X^, one can expect a considerable contribution from this mechanism in 
B —>• J/T-|-Xs due to the large ratio of the Wilson Coefficients between color-octet operator 
and color-singlet operator responsible for B —> J/T -|- X^. QCD sum rule0 so far is a 
powerful tool to deal with hadron matrix elements. Although there is an up to 30 percent 
uncertainty in this method, we can still give a useful estimation of the nonfactorization effect 
mentioned above. 

The effective weak interaction Hamiltonian at a scale /i is given by 

He,, = -^va,v;. Y. Ci(A)OiM, (1) 

V 2 
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in the conventional notation, O 2 = (sL 7 At^&L)(cL 7 ^^CL), Oi only differs from O 2 in the 
way of color indices contraction, and the former is referred as color-octet operator. On the 
factorization assumption, the matrix element of Oi dominates the decay B J/'^Xg while 
O 2 gives zero. 

The computation of the inclusive B J/'^Xg decay rate from the Oi contribution is 
performed by calculating the correlation function 


T = * jd‘xe-*^B(v)\T\jl(x)]Mmv))(-,r + ^), 

where = SL'jfj^bL. The corresponding decay rate is given by 

dr ^ 

PjiqidEjjqi 

where /j/^ is defined by ( 0 |c 7 ^c| J/T(p, e)) = ifj/^e^. 

At the leading order in the operator product expansion (OPE), ImT has the form 

2Ej/^?T7.ft 3Ej/qrMji\^ -|- 


( 2 ) 


( 3 ) 


ImT = 116 {ml + — 2mbEj/^)M^ 




( 4 ) 


It is well known that the next to leading order contribution in OPE is suppressed by 
Therefore, on the factorization approximation, (|^) almost gives the total 
decay rate of i? —J/TX^. The main non-factorization contribution arises from the interfer¬ 
ence terms of Oi and O 2 as shown in Fig.2, which is of order of AQ^^/m^[Q. The coefficient 
ratio of operator O 2 over Oi is C 2 /C 1 ~ 20. Therefore, it is reasonable to expect a large 
enhancement when contribution of Fig.2 is taken into account. 

The decay rate from Fig.2 is obtained as 

2Gl\V*Vtb\‘^Ci*C2,,_^^, dp^ 


dT = 


2M, 


B 




(2vr) 


( 5 ) 


3 








where and are defined as 

= I d^x^-^Ti5{p^ - Mj/^)(0|jlJ/vI/(p))(J/vI/(p)|j-(a;)jf (0)|0) 

= ImK^^i3o,{p) , 

rpvij.0a ^ / d^a;e“*^'^i(i?|6L(a;)7^S's(a;, 0)7 '^(—z(7sG^“)6l(0)|-B) , 

( 6 ) 

= 'xd{p^ - MJ|^){{]\J^\J/^!{p)){J/^!{p)\(^i^,{-igsG0a)c\^) 

= ImJ d'^xe^P'^i{0\j^{x)c{0)'yf,{-igsGpa)c{0)\) 

= dmKl^p^{p) , 

j.,ut.ga ^ Jd^xd^ye-^P-^t^{B\bL{y)YSs{y,x)rSs{x,0)rbLmB) 


In equation (^, we have assumed that other resonances contributions to I and /' are small 
compared with J/d' . Ss{x,y) is the propagator of s quark, = cq^c, j“ = cq^^c and 
ja5 _ ’^g only keep the leading order terms in OPE. Jm in the definition of / and 

I' should be understood as an operator which selects the imaginary part of the scalar form 
factors which factorize the followed matrix element. 

In this paper, / and /' are estimated by QCD sum rule. The matrix elements Klfi^a 

in equation are obtained by calculating the diagrams shown in Fig.3 and Fig.4, 


K, 




= + 


K' 


u^Pa 


OTXd 


(7) 


where is from Fig.3 and others are from Fig.4 and correspond to the diagrams with the 
indices a, b, c and d respectively. 

In a conventional way, the fixed point gauge x • y4(x) = 0 is taken. At the first order of 
the expansion in term of x^ is, 

Mk) = =^G.40)^'5-(A.). (8) 
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We obtain 


- ^ 4 y ^27,)D Qq 0 


Tr{ 


-la 


Iv 


k — rric k — q — me k — p — q — m, 


-757m + 


-7m ^ 


-757m^ 




Tr{ 


96 

-lo 


k — me k — q — me k + p — m, 

>/ ^ 


'li'}\g=0y 


{27i)^ dq^ dki dk^ 


-la 


k — me k — q — me k — q + k2 — m, 

ufifia 


■757m 


-lu 


k-b __*a A _A A \/ 

^ufifia ~ gg (<^pcr<5p/a-' 5p^i5ap'){ q gG ) ^Tr{ 2 2 2 2 


k + p — ki — me k — ki — me 

>/ " ^ ^ " 


1 p''\\q=ki=k2=0i 


(27r)'^ dqf^ dk^ dk^ 


Tr{ 


-7 


-7p' 


757m- 


-7a 


k — me k — q — me k — q + ki — me k + p — k2 — me k + p — m, 


'lu\\q=k\=k2=0i 


- V'' 5 v)( 97 ")irr{~ 77 } I 


Tr{ 


-lo 


k — me k — q — m. 

K^vpjla = -^{SerpSpe^ - 5„e.8p0){glG‘^) J 


d^k d d d 

{ 2 ti)^ dqf^ dki dk^ 
i i 

Iv 


-757m f , - r r la'j - r i i 7p'}|q=fci=fc2=0) 

; k + p — ki — k 2 — me k + p — ki — me k — ki — me 

d^k d i i i 

Tr{-f - ipj —:-7m r—:-}l9=o- 


{ 2 T^)Ddq^ ^k 


me k — q — me k — p — me 

(9) 


The computation of (§) is tedious, we will not display the procedure. Instead, we give the 
hnal result at the end of this paper. 

Using the dispersion relation 


Ki/pfia ip') 

---j 

f ( 


( 10 ) 


s — p^ 



^ v jj.ga^P') 

-L 

r 

S — p2 

'"bs. 

(11) 


and the fact that lup^a and are proportional to we get two sum rules 

lui^paip) = vr5(p^ - M^/^)Mj/^[^K^p0a{p)\p2^ol (12) 

^Ififiaip) = — Mj/^)Mj/^[-jj^Kl^p^{p)\p2^Q], (13) 


where = —p^. -^K^ppaip)\p' 2 =o must be understood that only derivatives of the scalar 
form factors of the matrix element are set at = 0. To choose = 0 is for the sake of 
convenience. In general, we can set p^ at any point in the range —p^ + m^ » ^^qcd 1° S®1 
a sum rule. In order to get rid of the dependence on the subtraction in the loop calculation, 
we have used the sum rules for derivatives of and instead of themselves. 
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The calculations of and are performed, in which we need the identities 


l^^{B{v)\bTgsGo,fih\B{v)) =-^Tr{T{l + v)a^fi{l + v)}, (14) 

(S(n)|6r6|5(n)) = ^Tr{(l + h)r(l + n)}, (15) 

where T is any kind of Dirac structure and /i^ is the value of the strength of the chromo- 
magnetic interaction of the b-quark inside in B hadron, 

= |(Afy - Ml) K aAGeV\ ( 16 ) 


After a tedious calculation, we arrive at 




- mbm^ - Sm^Ej/ii, AArribEj^^ - SSrribm^ - llm^Ej/q, (asC^). 


24:07r'^mf. 


12096m? 


247r 


ml + 61m^ — 6 mbEj/^ {agG"^) 


240m^ 


6 % 


(17) 


In the second equation of (pT]), we have already used partial integration. Using the standard 
numerical values 


(asG^) = 0.04G'eU^ m^ = l.dGeU, mb = A.bGeV 

Gi{mb) = (2C+(mb) - C'_(mb))/3 = 0.133, G2{mb) = G+{mb) + G_{mb) = 2.21, (18) 

fj/'ir = 0.38GeU, 

we obtain the nonfactorization contribution arisen from color-octet amplitude as large as 


5r(5 ^ J/T + X,) ^ 1.25ro(5 ^ J/T + X,), 


(19) 


where Tq is the result with factorization assumption. This result is sensitive to the value of 
me, if we choose m^ = 1.5GeU, the ratio would be ~ 0.7. The non-perturbative effect from 
the the gluon condense is much smaller compared with the perturbative diagram, so, the 
result is not sensitive to the parameter (cksG^). However, since the decay rate in (|) is inde¬ 
pendent of me, the non-factorization contribution is very large any way. One may notice that 
we have only taken account of the interference term shown in hg.2. If the non-factorization 
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effects from pure O 2 contribution are taken into account, which are very complicated to be 
calculated and will be presented in a separate paper |^, the enhancement will be hopeful to 
explain the discrepancy between the theoretical predictions and the experimental data. 
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Figure Captions 


Fig.l Nonperturbative effect in i? —+ 7 

Fig.2 Nonfactorization effect in 5 —+ J/'F 

Fig.3 The perturbative diagrams contribute to 

Fig.4 The gluon condense diagrams contribute to and 
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Fig.l 






